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Abstract 

We investigate the renormalization of "nonlocal" interactions in an effective field 
theory using dimensional regularization with minimal subtraction. In a scalar field 
theory, we write an integro-differential renormalization group equation for every pos- 
sible class of graph at one loop order. 
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1 Introduction 



In its traditional form, an effective field theory calculation goes like this: Start at a very large scale, 
that is with the renormalization scale, fi, very large. In a strongly interacting theory or a theory with 
unknown physics at high energy, this starting scale should be sufficiently large that nonrenormaliz- 
able interactions produced at higher scales are too small to be relevant. In a renormalizable, weakly 
interacting theory, one starts at a scale above the masses of all the particles, where the effective 
theory is given simply by the renormalizable theory, with no nonrenormalizable terms. The theory 
is then evolved down to lower scales. As long as no particle masses are encountered, this evolution 
is described by the renormalization group. However, when \i goes below the mass, A, of one of the 
particles in the theory, we must change the effective theory to a new theory without that particle. 
In the process, the parameters of the theory change, and new, nonrenormalizable interactions may 
be introduced. Both the changes in existing parameters, and the coefficients of the new interactions 
are computed by "matching" the physics just below the boundary in the two theories. It is this 
process that determines the relative sizes of the nonrenormalizable terms associated with the heavy 
particles. 

Because matching is done for /i m A, the rule for the size of the coefficients of the new operators 
is simple for fi « A. At this scale, all the new contributions scale with A to the appropriate power 
(set by dimensional analysis) up to factors of coupling constants, group theory or counting factors 
and loop factors (of 16ir 2 , etc.) ]!]]. Then when the new effective theory is evolved down to smaller fj,, 
the renormalization group introduces additional factors into the coefficients. Thus a heavy particle 
mass appears in the parameters of an effective field theory in two ways. There is power dependence 
on the mass that arises from matching conditions. There is also logarithmic dependence that arises 
from the renormalization group. 

The matching correction at tree level is simply a difference between a calculation in the full 
theory and a calculation in the low energy effective theory 

/ 8£°{9) = S° Ch+c ($) - S° c {9) 

J (1) 

If virtual heavy \ ( rf\\ 

I \ particle trees J \ / 

where +c ($) denotes the light particle effective action in the full theory and S^(^) denotes 
the same in the low energy theory 0. The matching correction so obtained is nonlocal because 
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it depends on p/A through the virtual heavy particle propagators. It is also analytic in p/A 
in the region relevant to the low energy theory, i.e. for characteristic momentum << A. Thus 
it can be expanded in powers of p/A with the higher order terms decreasing in importance: this 
corresponds to a local operator product expansion in the domain of analyticity, equivalent to a local 
nonrenormalizable Lagrangian which can be treated as an honest-to-goodness local field theory. 
However in general an infinite series of terms of increasingly higher dimension are generated by 
matching at tree level. These cause no problem when the scales are well separated, because their 
effects quickly become negligible. But if there are two or more scales close together (such as m t and 
Mz may be), then we may not be justified in ignoring terms at higher orders in the expansion. How 
does one understand how to interpolate smoothly between the well-understood situation in which 
the scales are very different and the well-understood situation in which the scales are very close 
together? How can we keep track of all the infinite number of higher derivative operators efficiently? 
Is it possible to deal with the nonlocal effective Lagrangian directly, without expanding? In the 
particular context of a scalar field theory, we will attempt to answer questions at the one loop 
level. The approach will be direct. We will manipulate the nonlocal interactions as if they are 
expanded in a momentum expansion, and then show that the resulting f3 functions for the terms in 
the momentum expansion can be collected into integro-differential renormalization group equations 
for the nonlocal couplings. 

The present paper is organized as follows: in section 2 we calculate the (3 function for a nonlocal 
four point coupling arising from a one-loop graph with two internal lines in nonlocal (ft 4 theory. 
We give this example before the general results of the following sections in order to point out the 
important features of our method. In section 3 the method is then applied to obtain the contribution 
to the renormalization group equation for a general class of graphs in an arbitrary massive nonlocal, 
non-renormalizable parity invariant scalar field theory. 

2 Basic Example 

The matrix elements in a general nonlocal field theory || are calculated by writing down a La- 
grangian with 'smeared' vertices. For instance, for a nonlocal (ft 4 interaction in a scalar field theory 
with the discrete symmetry (ft — > —(ft, the interaction term in the action is 

5*4 = J dxidx 2 dx 3 dx4F (xi — x 4 , x 2 — x 4 , x 3 — x i )(ft{xi)(ft{x2)(ft{x' i )(ft{xi) (2) 
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Figure 1: Basic nonlocal 4 vertex. 

where Xi are spacetime coordinates and F is a nonlocal 'form-factor'. Energy- momentum conser- 
vation at each vertex of the corresponding Feynman graph is expressed in terms of the Fourier 
transform G of F: G = G(pi,P2,P3) for the basic 4 interaction shown in figure 1. Bose symmetry 
implies that the nonlocal coupling G is symmetric and satisfies 

G(pi,P2,Ps) = G(p!,p2,-pi -P2-P*)- (3) 

Of course, Lorentz invariance dictates that in the final expression for the matrix element in mo- 
mentum space, only scalar products of momenta will appear as arguments of G. Implicit in the 
definition of G is a mass scale A which sets the limit for the region of analyticity of G, and for 
characteristic momenta p < A, G is analytic. We go to the effective low energy theory by expanding 
in a Taylor expansion in p/A. One can think of this Taylor expansion as the formal implementation 
of a local operator product expansion of G. At this point, in going to the effective low energy 
theory, we have actually changed the high energy behavior of the theory so that integrals which 
were convergent in the full theory are divergent in the effective theory. This trades logs of A in the 
full theory for anomalous dimensions in the low energy theory. This trade allows us to calculate 
the logs more simply and to sum them using the the renormalization group. 

Let us first discuss as an example the renormalization of the nonlocal 4 interaction from the 
graph shown in figure 2. For massless fields, this is the only contribution in one loop. In the 
massive case, there is also a contribution from a tadpole graph. In addition, in either case, there 
are renormalization of <f) 2k interactions for k > 2. We will systematically consider them in the next 
section. 
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Figure 2: Feynman graph contributing to the renormalization of G. 

The Feynman integral for the graph in figure 2 is 

d 4 k G(p 1 ,p 2 ,k)G(p 3 ,p±,k+pi+p 2 ) 



Combine denominators 



(27r) 4 [k 2 — m 2 + ie] [(k + pi + P2) 2 — rn? + ie] 



1 

[k 2 — m 2 + ie] [{k + p± + P2) 2 — rn 2 + ie] 
da 



(4) 



(5) 



and shift momenta 



to obtain 



[(1 — a)k 2 + a(k + p\ + P2) 2 — m 2 + ie] 2 

k = £- a (p 1 +p 2 ) (6) 



Now, here is the crucial point. We get the low energy effective theory by expanding 
the Gs in a momentum expansion. An equivalent procedure is to treat the Gs as if they 
were analytic everywhere in momentum space. Thus we can deal with the I dependence of 
G by writing a symbolic Taylor expansionS Here we are effectively doing the momentum 



1 Analyticity of the nonlocal couplings may alternatively be exploited by formally Laplace transforming in the 
loop dependent arguments, which yields instead of e l '^G an exponential e~ e ' s G in terms of the Laplace transform 
of G, along with an extra Laplace inversion integral. The rest of the computation is essentially identical. 
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expansion. But the key is that we can resum the final result into an finite integral over 
the original nonlocal couplings. 



G(pi,P2,£- a(pi +p 2 )) G(p 3 ,p 4 ,£+ (1 - a)(pi + p 2 )) 

(8) 

= e e * [G( Pl ,p 2 , q) G(p 3 ,p 4 , q + pi+ p 2 )] q= _ a{pi+p2) 

Now all the dependence on the loop momentum is in the denominators and in the exponential, 
so we just have to do the Feynman integral of the exponential 

f*L (9) 

J (2tt) 4 lp + a (l-a)(p 1 +p 2 ) 2 -m 2 + ie\ 2 

or to be more precise, the dimensionally regularized integral 

d A ~ e £ 

(27r) 4 - e /i- 2e \p + a {l - a)( Pl + p 2 ) 2 - m 2 + ief ( ' 

We will do this by manipulating the exponential like a power series, because the basic assumption 

is analyticity in the momenta. It is at this point that we have irrevocably changed the high energy 

behavior and gone over the effective low energy theory. Because of symmetry, only even terms in t 

contribute, 



f i r (^r f:l:n 

r tr {2r)\ J (271)*-^-^ [i2 + a (l-a)(p 1 +p 2 ) 2 -m 2 + ie] 2 



We calculate this as follows 



^£ £ ■■■£ 



(2r)!/(2 r r!) terms (12) 



where 



W2 ■ ■ ■ K-y* + P erms 



n = J dQ e . (13) 

Contracting with g^ 1 ^ 2 gives 

A r _! = (4 + 2(r-l)) A r = 2(r + l)A r (14) 

or 

A r = —, rr (15) 
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and thus 



So the integral is 



1 r dVL f / d 



2r 



(2r)! J \ dq 



1 




r 




4 r r!(r + 1)! 









(16) 



E 

r=0 



1 



d 



4-ei 



I? 



2l r 



4 r r!(r + l)! J (2ir)^^ ^2 + a ^ _ a ^ pi + p ^2 _ m 2 + 



(17) 



which we can calculate in the usual way. Wick rotate 

(-1Y r d A ~H 



00 

*E 



f A r r\(r + 1)! J (2tt) 4 -^-^ [£ 2 + A(a) 2 ] 2 
where 

A(a) 2 = -a(l - a)(pi + p 2 ) 2 + m 2 , 

and A(a) 2 is positive for Euclidean momenta. Now doing the i integral gives 

2" 



2i 



(4vr) 2 e/ u- 2e 



4 r r! 2 



2r 



r=0 



A, 



+ 



(18) 



(19) 



(20) 



where we have dropped everything but the 1/e pole and the associated \i dependence. We can write 
the result as 



2i 



d 



{An) 2 eiJ,~ 2e dx 



E 



X 



-xA(a 



,2lr 



A, 



(21) 



x=l 



^ 4 r r!(r + 1)! 

This can then be turned back into the (Euclidean) angular average of an exponential, just inverting 

(HBD, 

(22) 



(47r) 2 e/i~ 2e cte 



x e 



x=l 



where e is a Euclidean unit vector. Now that we have done the £ integral and extracted the 1/e 
pole, we can use the Taylor series to put the form (|22"D in terms of the Gs. The result is 



d 



(87r) 2 e / u~ 2e dx 



dn e /-i 

-— r- / r/n ./• 



G(pi,p 2 , y/xA(a)e - a(p 1 + p 2 )) 



G(p3, ~Pi -P2 ~P3, VxA(a)e + (1 - cv)(pi +p 2 )) 



(23) 



x=l 



with the derivative evaluated at x — 1. Thus we have reduced the divergent part of the original 
Feynman graph to an integral of the Gs and their derivatives over finite ranges (for Euclidean 
momenta). These last integrals are well-behaved and could be done numerically. Hence, with 



Figure 3: Feynman graph contributing to the running of 4> 6 coupling. 



Q = 2ir 2 in for Euclidean four space, we obtain the one-loop (3 function for the nonlocal interaction, 
G: 

1 d 



dVt P I dax 
o 



167r 4 dx 

G(p 3 ,-Pi -p 2 -P3,VxA( a ) e + (1 -a)(pi +P2)) 



G(p 1 ,p 2 , y/xA(a)e - a(pi + p 2 )) 



(24) 



+ 



cross 



x=l 



terms 



Now in the massless case, where there are no other contributions to the renormalization of the 
4 coupling, the "running" nonlocal coupling satisfies the integrodifferential renormalization group 
equation 

ft 

(25) 



A useful check on this result is obtained by going to the 'local limit', i.e. G — » g where g is the 
usual coupling constant of local 4 theory. Indeed, we know that the contribution to the (3 function 
in local (f) A theory at one loop is 

(26) 



Hg) ^ 0. 



16^ 

Substituting G = g, inserting an explicit symmetry factor of 1/2, and summing over the crossed 
graphs, (PH) indeed yields 

( 3 G(p 1 , P 2,P3) as G -> g. (27) 

Now the nonlocal quartic coupling G can induce changes in the coupling terms with more fields. 
For instance, the one-loop diagram relevant to the renormalization of the </> 6 coupling is shown in 
figure 3. 
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Figure 4: Feynman graph that vanishes for massless (j)s. 



For massless fields, this does not mix back into G in a mass independent renormalization scheme, 
because the relevant Feynman graph of figure 4 vanishes. However, in the massive case, the graph 
of figure 4 gives a non-vanishing contribution. We will now systematically compute the /3 functions 
of the general nonlocal couplings in a massive nonrenormalizable theory. 

3 General Results 

Working within the formalism of a massive nonlocal effective theory induced by some unknown 
full theory, we are forced to consider an effective Lagrangian with operators with an arbitrary 
number of low energy fields. In this section, we will compute the one-loop running of a general non- 
renormalizable, nonlocal scalar effective theory with <fi — > — symmetry. Specifically, our purpose 
is to isolate the dimensional regularization pole of a 2m point function of type n (i.e. with n 
vertices or propagators) at one loop. The plan is as follows: we first construct the expressions 
corresponding to the Feynman integral of a general graph with a special choice of momentum 
labelling conventions. We then describe the application of the method of the last section to two 
relevant examples: the tadpole renormalization of the nonlocal '0 4 ' coupling with a 6 operator 
insertion, and the renormalization of a nonlocal 6 coupling with three 4 operator insertions. The 
last example is useful in highlighting the special features of renormalizing graphs with more than two 
internal lines. Finally, we attack the general case, and give the complete expression for the 1/e pole 
of the 2m point function of type n at one loop as a surprisingly compact multi-Feynman-parameter, 
multi-dimensional angular integral. 
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Figure 5: Type-n Feynman graph. The blobs signify arbitrary 2r insertions. 
3.1 Preliminaries 

We assume that the dimensionally continued d = 4 — e dimensional non-local Lagrangian has a 
— > —(f) symmetry and hence has an interaction term proportional to 

oo 

Ant. =E^ (r " 1)G 2,0 2r . (28) 
r=l 

The mass scale fi is introduced to keep the dimensions of the nonlocal couplings G<i T fixed under 
dimensional continuation. The interaction is not normal ordered!, and each G^r is some nonlocal 
function which is analytic in the region under consideration, depends as a consequence of momentum 
conservation on 2r — 1 linearly independent momenta, and may have dimensions proportional to 
some power of an implicit scale of nonlocality A. 

A diagram with n internal lines will be called type-n. At one loop, a type-n graph has n vertices 
connected to external lines. Let 

n 

N = ^Tvi, (29) 
i=i 

where Vi is the number of external lines emanating from the ith vertex. The loop integral for the 
type-n one-loop renormalization of the A" point function (shown in figure 5) is 

d*k Ui=i Gi i+2 {- --,k + Qi-i) 



(2TiY^)(k + Q l f...(k + Q n ^f 
where the upper index on the G's is used to number the vertices as one goes along the loop. 



2 Hence fields at the same point can be contracted, and tadpoles will occur explicitly. 
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Here all external momenta are taken to be incoming (signified by . . ., in the numerator, different 
for different G l ) and we have defined 

Qo = (31) 

Qi = Pi 
Q2 = P2 + Q1 
Qs = P3 + Q2 



j 



Qi — Pi + Qi-i — /.Pi 

Qn Qn—1 "i" Pfi 0) 

and Pi is the sum of the external momenta flowing into the ith vertex. Energy momentum conser- 
vation is simply Q n = 0. Symmetry factors are suppressed. 

Using this notation, the loop integral for the for the 4 coupling that we computed in the last 
section is obtained by setting n — 2, vi — 2, v 2 — 2: 

r d 4 k Gl{p 1 ,p 2 ,k)Gl{p 3 ,p A ,k + Q 1 ) 

i(27r) 4 (k 2 + ie)[(k + Q 1 ) 2 + ie] ' 1 } 

Using the generalized Feynman identity 

1 r(pi + p 2 + • • • + pn, 



A?A p 2 2 ...A p n T( Pl )T(p 2 )...T(p n ) 

a pi _1 o;f _1 . . . <" -1 <5(l - ai - a 2 - . . . - a n ) 



(33) 



x / doii . . . da T 
Jo 



(aiAi + a 2 A 2 + . . . + a n A n )p^+-+P- 
combine the denominators 



(k*)(k + Q 1 ) 2 ...(k + Q n - 1 y 



(34) 



and shift the loop momentum 

71—1 

fc = * - £ Qs®s (35) 



s=l 

so that 



n— 1 n— 1 

^ 2 = ^ + (E ^«*) 2 ( 36 ) 

s=l s=l 
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which, when put into the denominator gives 



n— 1 n— 1 n— 1 

D = [£ 2 + (^Q s a s ) 2 -2£-£Q s a s -£ 2 5> t (37) 

8=1 S=l t=l 

71— 1 71—1 71—1 77,-1 77,-1 

-(E Q,« s ) 2 £ a* + (2* • £ Q-a.)(£ a*) + E 

s=l t=l s=l t=l r=l 

71—1 71—1 71—1 71—1 71—1 



+ E (E QsUs) 2 oi T - 2 E * • E &«»«r + E Q 

r=l s=l r'=l s=l r= 

n— 1 Ti—1 Ti—l 

+2 £ £ • Q r a r - 2 E Qr(E Q s a s )a r } n . 



r=l r=l s=l 

After cancellation ( |37"D yields the denominator 

D = [£ 2 — A 2 ] n (38) 

where 

71—1 71—1 

A 2 = (J2Qsa s ) 2 -J2Qrar- (39) 



s=l r=l 

For the massive case, we simply substitute in fl53| ) 

A; 2 -> A; 2 -m 2 (40) 
(A; + g r ) 2 -> (A; + g r ) 2 -m 2 (41) 

so that we just get an additional term in ( [T7|) when combining denominators after the shift (|35|): 

71— 1 71—1 

- m 2 (l - E a s ) + £(-m 2 K = -^ 2 - (42) 

8=1 S=l 

So the general form for the denominator of the right hand side of ([34]) with massive fields is 

71—1 71—1 

D = [k 2 — (E Q s a s ) 2 + (E Q> r ) - m 2 f. (43) 

s=l r=l 

Now the shift fl35|) changes the argument of the numerator factors in ( |30"D also, to yield the final 
expression for the Feynman integral! 



3 The reader may check that with n — 2, «i = V2 = 2, the above relations give the correct integral for the example 
of the first section: 



"" ' ( 2 ^) d [£ 2 + a(l-a)Q2_ TO 2] 2 



with Qi = pi +P2- 
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Important remarks on notation : (1) In intermediate steps of the computations, we denote the 
dependence of each nonlocal function on the (distinct!) external momenta by ellipsis. This is useful 
since the external momenta play a trivial part in the loop integration, and may be reinstated by 
examination in the final expressions. (2) Also, since at one loop any vertex shares two and only two 
lines with the loop, by energy momentum conservation the loop momentum appears only once as 
an argument of any G l and will be put in its last slot. (3) The first vertex will have only the loop 
momentum as its last entry. 

Crossing : The external momenta can be exchanged amongst themselves. The final result, 
however, depends only on the number of distinct Lorentz invariants of the four momenta Pi (i = 2m 
for renormalization of the 2m point function) under the condition Yl1=\Pi — an d pj = m 2 . This 
equals the total number of graphs related by 'crossing' where the crossed graphs can be obtained 
by exchanging external momenta. We will give explicit expressions only for one member of each 
crossed set. 

Counting i's : Ignoring for the moment the i's appearing due to Wick rotation and integration 
(see below), the integrand itself gives no powers of i. This is seen as follows: each propagator is 
2_^2 +ie , and each nonlocal vertex has a Feynman rule —iG<i r - Since the number of vertices equals 
the number of internal lines at one loop, for a type-n graph we obtain (—i) n i n = (—l) n (i) 2n = 
(— l) 2 ™ = 1. So the only source of i's and minus signs is the integration formula. 

Computational Tools : 

In order to do the dimensionally regularized loop integrals we will need the well-known integra- 
tion formula (for Euclidean momenta): 



/ 



d «tf- A *)»- m { l) {A) r(f)r(n) ' (46) 

and the expansion 

r(-n + e) = ^[l + (l + U... + i- 7 )+(9(e)], (47) 

where 7 = 0.5772157 is Euler's constant. 

Also note that in 2k — 2 dimensions (k a positive integer) 

" .£ £ ...£ 



(2r)!/(2 r r!) terms, (48) 



2k-2 



W2 ' • • 9nar-lVOr + P erII1S 
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where we have defined 

Contracting with g^^ 2 gives 

A 2 ^ 2 = [2k-2 + 2{r - 1)] A 2k ~ 2 = (2k + 2r - 4) A 2 r K ~ 2 (50) 

So 

A 2k ~ 2 = — (51) 

and thus in 2k — 2 dimensions 



1 dtt 2 e K - 2 




, d 


\ 2r 


(2r)\ n 2 - 2 


(< 


dq) 


1 




r 


{dq) 


4 r r!(r + k - 2)\ 









(52) 



The importance of the result of (|52|) will be seen when we isolate the pole pieces appearing from 
graphs with more than two internal lines. The point is that the leading terms in the Taylor 
expansion then give Feynman integrals which are manifestly convergent by power counting, and 
divergences appear only at some higher order in the Taylor expansion. This implies that the sum of 
dimensional regularization poles does not begin at zero, and we cannot immediately write the result 
as an angular integral over a Euclidean unit vector in four dimensions, as we did in the last section. 
One choice is to redefine the infinite sum to start at zero and compensate for the redefinition by 
subtracting off a finite sum (which vanishes under the action of differential operators in the dummy 
parameters). Alternatively, using (0), we can re-sum the poles for a graph with n > 2 internal 
lines, without reference to any dummy variables , in terms of a 2n — 2 dimensional angular 
integral! The example of the 2 — 2 — 2 graph below will give the explicit details of how this is done. 

The 6 operator can be renormalized at one loop by a type-1 (tadpole) with a Gg coupling 
insertion, a type-2 graph with Gq^G^ insertions, or the most convergent graph, of type-3, with 
three G4 insertions. For the last one, the relevant graph with incoming external momenta is given 
in figure 6. 

As mentioned in the last section, for massless fields this will not mix back into G4 in a mass 
independent renormalization scheme, because the Feynman graph in figure 7 vanishes. However, 
for massive light fields the tadpole graph does not vanish, and we will therefore evaluate this type-1 
'tadpole' graph first. 
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Pi P2 P3 Pa 
Figure 7: Tadpole graph that vanishes for massless 0s. 
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3.2 Tadpole diagram contribution to /3q 4 

We need to do the integral (see figure 7), where r = 1 in (|28|) , and n = l,Vi = 4,p 4 = — (P1+P2+P3 

2€ ! d d k G e (p u p 2 ,P3,P4,k) 



J (2n) d (k 2 -m 2 + ie) ' 1 ' 

No momentum shift is needed since the denominator is purely quadratic in the loop momentum. 
Now because Gq is analytic, we can deal with the k dependence by writing a symbolic Taylor 
expansion 

Gs(pi,P2,P3,P4,k) = e k ^ [G 6 (p 1 ,p 2 ,p 3 ,p4,g)]g = o (54) 
Then we just have to do the dimensionally regularized integral 

d*- € k 



1 



e k 8q 



(27r) 4 - e /i- 2e [k 2 -m 2 + ie] 



(55) 



We again do this by manipulating the exponential like a power series, because the basic assump- 
tion is analyticity in the momenta. Because of symmetry, only even terms in k contribute, so the 
integral to be done equals 



E 



1 



d^- e k 



J Q (2r)! J (2tt) 4 ~^~ 2 " [k 2 -m 2 + ie] 
which finally yields the integral (just as in example of last section) 

1 r d^k \- k 



(56) 



21 r 



E 



9.V 

dq) 



(57) 



o 4 r r!(r + 1)! J (27r) 4 -^- 2e [k 2 - m 2 + ie] ' 
After performing a Wick rotation, and doing the integral using ( (4"6| ) and ( f4T|) we obtain for the 1/e 
pole and the associated \i dependence: 



im 



E 



1 



\m 



2ir 



+ 



(5c 



(47r) 2 e/i" 2e ^0 4 r r!(r + 1)! 
This can then be turned back into the (Euclidean) angular average of an exponential, by just 
inverting flTB]) , 

im z r dilt m ,JL 

(59) 



me-i- 

6 d i 



(47r) 2 e/i- 2e J ft 4 

where e is a Euclidean unit vector, and f2 4 = 27r 2 . Now that we have done the loop integral and 
extracted the 1/e pole, we can use the Taylor series to put the form ( 59" ) in terms of Gq. The result 
for the pole part then equals 



im 



32vr 4 e/i 



-2c 



J dQj [Ge(pi,P2,Ps, -Pi -P2 -P3,me))] 



(60) 
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Thus we have reduced the divergent part of the original Feynman graph to an integral of Gq over 
a finite region (for Euclidean momenta). Hence we have obtained the tadpole contribution to the 
one- loop j3 function for the nonlocal interaction, G^: 



^G m (pi, P2 ,ps) = ^-J^G 4 [4](p 1 ,P2,P3) (61) 



2 

m 



J dtt*GQ(p 1 ,p 2 ,p 3 ,-p 1 -p 2 -p 3 ,me), 



32vr 4 

where the subscript on G on the LHS denotes the coupling that is renormalized, with the number 
of external lines at each vertex for the graph under consideration in square brackets (i.e. type of 
graph considered) .@ 

3.3 2-2-2 diagram contribution to (3g 6 

The Feynman integral from figure 6, with Qi = p\ + P2, Q2 = Pi + P2 + P3 + Pa, J2iZiPi = Oj an d 
V\ = v 2 = v 3 = 2, equals 

3e f d * k G4(p 1 ,p2,k)G 4 (p 3 ,p i ,k + Q 1 )G4(p 5 ,p 6 ,k + Q 2 ) , , 

_/i J (2vr) 4 (P - m? + ie)[(k + Q1) 2 - m? + ie][(k + Q 2 ) 2 - m? + it}' 1 ' 

Combining denominators and making the shift k = i — [Q\Oi\ + Q 2 a 2 ), we get the denominator 



D 



- Q^i^ai - 1) - Q 2 a 2 (a 2 - 1) - 2Qi • Q 2 a\a 2 - m (63) 
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so that the integral becomes (with a factor of 2 from the Feynman trick) 

= 2jY[d ai J -^J- i ^[G i (p l ,p 2 ,£-(Q l a 1 + Q 2 a 2 )) (64) 
Ga(P3,PaJ- (<5i(«i - 1) + Q2tt2))G , 4 (p 5 ,p 6 ,£-Qiai - Q 2 (a 2 - 1)) 
Define the numerator 

N(£,p,a) = [G A {p x ,p 2 ,t- (Qiai + Q 2 a 2 )) (65) 
G 4 (P3,P4,^ - (<5i(ai - 1) + Q 2 a 2 )) 
Ga(P5,P6,£- Qioci - Q2(a 2 - 1))]. 



4 Also note than in the local limit, replacing Ge by a 'constant' ge, we get /3 94[4] = " 6 ^| which may be directly 
obtained in the local limit as a check. 
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Now, because the Gs are analytic, we can deal with the I dependence of G by writing a symbolic 
Taylor expansion 



N{£,p,a) 



e d i 



G 4 (pi,P2,q) 
C 4 (P3,P4,g + Qi) 



(66) 



-(Qiot 1 +Q 2 a 2 ) 



and we just have to do the dimensionally regularized integral 



e 9 i 



(27r) 4 - e /i- 3e [P -A 2 + ief 
where A 2 = Qja^^ - 1) + - 1) + 2Q 1 ■ Q 2 aia 2 + m 2 . 



(67) 



Again, we do this by manipulating the exponential like a power series, because the basic as- 
sumption is analyticity in the momenta. Because of symmetry, only even terms in I contribute, and 
doing the t integral yields a sum for the pole piecesi: 

°° 1 r / ^ \ 2_| 

E t^ttt ^, (^ r - J 



d_ 



(68) 



8ttV~ 3€ r =i4 r r!(r-l)! 

where we have retained only the 1/e pole and the associated \i dependence. We cannot yet convert 
this sum to an angular integral, because it starts at r = 1, which just reflects the fact that the 
leading term in the Taylor expansion gives a convergent Feynman integral, with no 1/e pole piece. 
To put the contribution of the sum back into the G's, we attempt to massage it further: 

2" 



E 



1 



{A r r\(r- 1)! 



(A 



2\r-\ 



d_ 



dx 2 

dx 2 

cP_ 
dx 2 



E 



\ 4 r r!(r + 1)! 
1 



2\r-l x r+l 



X 



J 4 r r!(r + 1)!M 2 



Y- 

jL^i At 



X 



{^)^) r x r 



dq) 
[dq, 



x=l 



X 

A 2 



2 = 1 



(69) 



x=l 



n A r r\(r + l)VA 2 ' 

Note that the finite sum subtracted off to redefine the sum to start at zero gets annihilated by the 
differential operator. 



3 Note also that the factor of 2/e gets compensated by a factor of 1/2 from a T function 
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The terms in the last sum (B^) are exactly of the form seen before in flTB]), and it can be written 
in terms of a four dimensional Euclidean angular integral over a finite range: 



dx 2 



x 



A 2 J ft 4 



(70) 



x=l 



which can be put back into the G's by inverting the Taylor expansion. The contribution to the f3 
function for the nonlocal 6 vertex arising from a graph with three G4 vertices is thus 



0, 



G 6[2,2,2](Pl>— Ps) 

1 d 2 



(71) 



WJ daida2 A 2 



x=l 



8tt 2 dx 2 

Ga(pi,P2, ^/xAe - (Qiai + Q 2 «2)) 
Gi(p 3 ,P4, VxAe - {Qxotx + Q 2 a 2 ) + Qi) 

Gi(p 5 ,Pe, \/xAe - <5i(«i - 1) - Q2(a 2 - 1)) 

cross 
terms 

where e is a Euclidean unit vector and the square root of A 2 is real for Euclidean momenta, so we 
have written it as A. 

Alternatively, we can eliminate reference to the extraneous parameter x by the following trick. 



Consider the sum of (p8|) again: 



E 



2\r-l 



{ 4 r r!(r - 1)! 



(A 2 y 



Pi, 



With 



we obtain 



p = r — 1 



E 



2\P 



9q j 



P+i 



4 [d q/ 
4 [dq, 



E 



d^t Ae-S- 

e 9 i 



(A 2 y 



dq) 



(72) 
(73) 

(74) 
(75) 
(76) 
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and using (|49| ) for four dimensions, / dQ A = 2tt 2 , we get a more compact form for the nonlocal (5 
function: 



(77) 



-^-j y dVt A e J da\da 2 Gi(pi,p 2 , Ae + q) 



64tt 4 \dq 
G 4 (p 3 ,p4, Ae + g + Qi) 

G , 4G°5,P6,-4e + <? + (22) 
cross 
terms ' 



-(Qiai+Q 2 a2) 



The last form is free of extraneous parameters. 

We reiterate that the evaluation of the 2 — 2 — 2 diagram above highlights two important features 
which will be crucial for the application of the method to the general case below. First, for graphs 
with more than two internal lines, the integral of the leading (zeroth) term in the formal Taylor 
expansion does not give a divergence (it is power-counting convergent), so the infinite sum for the 
poles starts only at some higher order: for a type n graph (n > 1) the first divergent contribution 
comes from the (n — 2)th term in the Taylor expansion, so the sum of poles starts at n — 2. To 
rewrite this sum in terms of the angular integral over a Euclidean unit vector, we can subtract off a 
finite sum which vanishes under the action of the differential operator of order n — 1 in the dummy 
parameter x. In fact, we can do much better - the parameter x can be eliminated by writing the 



pole sum a power of y§^j acting on a 2n — 2 dimensional Euclidean unit angular integral of 
G"s. Let us explain why this makes sense: For one-loop graphs with one or two internal lines, the 
four dimensional Feynman integral has a divergence which begins with the zeroth order term in the 
Taylor expansion (and higher terms in the Taylor expansion are even more divergent), so the sum 
of poles starts at zero. However, for graphs with more than two internal lines, enough powers of 
momentum in the numerator are required to cancel the denominator powers - thus the pole sum 
starts at a higher order. This sum can be made to start at the zeroth order by defining the 
loop integral in the appropriate number of higher dimensions. Hence the resummation of 
the pole terms for a sum starting at zero gives a higher dimensional angular integral. 
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3.4 2m-point function of type n at one loop 

For the non-trivial general graph with n > 2, after combining denominators and shifting the loop 
momenta, we geti the integral 

f l n-l „ idf i n n-l 

(n-1)! / Udaj l^ dT) \{^ {l) Gl i+2 {... ) i-Y,Qs^ + Q l -i) (78) 

J ° j=l Jl \ Z7T ) U i=\ s=l 

where Vi is the number of external lines emanating from the ith vertex, ctj are the Feynman param- 
eters, 

D=[f- CE QsVs? + (E Q>r) - m 2 T (79) 

s=l r=l 

is the combined denominator, and 

A(i) = e - Vi (80) 

carries the renormalization scale dependence. The ellipsis denote the dependence of each G l on the 
external momenta and the momenta Qi are defined as in (0)1 

For the renormalization of a 2m point function with n internal lines 

i=n 

$>i = 2m ( 82 ) 
i=i 

so summing up the fi dependence from all the vertices gives the overall power of fi appearing as 
fi Am where 

i=n ^ 

A m = J2vi(-)=me. (83) 
i=i 1 

Now, using again the fact that all the G's are analytic, we can write 

n n— 1 

l[G i Vi+2 (...,l-J2Qs<Xs + Qi-i) (84) 

1=1 s = l 



e d i 



nGt i+2 (...,?+Qi-o 

.i=i 



<?=-ELi 1( 3 s " 



3 The overall factor of (n — 1)! comes from the Feynman trick. 



7 By momentum conservation at the last vertex, we may make the structure of the last vertex simpler by making 
the replacement 

n — 1 n— 1 

G "„+ 2 (- + Qn-i) -> c: n+2 (...,-i + Yl Q><*>)- ( 81 ) 
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Now we just have to do the dimensionally regularized integral 



e-2- 

e a i 



(2tt) 4 -^- a ^ \t 2 - A 2 + ie] 



where 



n-l 



n-l 



a 2 = (EW-(E^) + 



m 



15) 



(86) 



s=l 



r=l 



Remembering the analyticity in momenta, we do this as usual by manipulating the exponential 
like a power series, with only even terms in I contributing. Doing the integral as before using the 
key formula (pEB]), and expanding the T function using (f4?j), with 



T(n - r - 2 + -) 
v 1' 



_(_l)2+r-n( 2 + r _ n y 

e 



(87) 



for n < r + 2, we get the pole piece 



(n-l)\- 



2,,-A„ 



57T /I 



E 

r=n— 2 



r^y+r-n(^yi+rf£2\2+r-n 

(n - l)!4 r r!(2 + r -n)! 



Quite nicely, the (n — 1)! coming from the integration formula cancels with the (n — 1)! appearing 
from the Feynman trick. The sum (|8~8"D can then be written as 

2" 



E 



(A 



2\2+r-n 



8vrV" Am e r ^l 2 4rr! ( 2 +r - n)! 



<9 



n— 1 



87T 2 /i~ Am e <9x n_1 

8ir 2 fi~ Am e dx 71-1 



oo 

E 

r=n-2 



(A 



4 r r!(r + 1)! 

(A 2 ) r x r 



'A' 
A, 



x=l 



E 

r=0 



2^n-2 ^ 4r . r ,/ r + ^ 



T 



(89) 



where we have defined the finite sum 

T 



x 



\n-2 



n—3 

E 



(A 2 ) r x r 



4 r r!(r + 1)! 



Pi, 



(90) 



(A 2 ) r=0 

taken to be vanishing for n < 3. As in the last example, this finite sum is subtracted off to make 
the first sum start at r = 0. Now the differential operator is of order n — 1, and the finite sum is of 
order n — 2, so the finite sum gets annihilated by the differential operator, giving finally the pole 
piece: 

x ^ (A 2 ) r x r 
a 

r=0 



i d 11 ' 1 



(A 2 



n-2 E A r r U r + !)! 



(91) 
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Now, inverting as usual using (|T6| ) for a Euclidean unit vector, we may cast this pole contribution 
as an integral over a finite four-dimensional Euclidean angular region: 

i d n ~ l 



Pole 



8vr 2 e/i~ me dx" 1 ' 1 
[ f 1 j TT 1 

y !^ y y j (a 2 )™- 2 



(92) 



x 



71— 1 



IJ G;. +2 . . . , y^Ae - Q s a s ) + Qz-i 

i=l V s=l / 

which gives the ft function in the formi 



J x=l 



ft 



G2m[t) 1 ,..., t)n ](Plv>P2m-l) 



(94) 



1 d 



n-l 



16vr 4 dx n - x 

i=n / 



1 j=n-l 



7=7i 

H 

T Vi+2 
i=l 

cross 
terms' 



71-1 



A/xAe - <5s«s) + Qi-i 



8=1 



X=l 



However, we can again do much better, and obtain a more compact expression in terms of a higher 
dimensional angular integral. Referring back to flS8|), we can write 

2" 



E 



in terms of 



for 



r=n _ 2 4T!(2 + r-ra)! 

p = r — n + 2, 
n > 3, 



d_ 



(95) 

(96) 
(97) 



8 With the property that at the last vertex in the product we may make the replacement, by momentum conser- 
vation, 

n-l 



G "„+2(- • • ) \fxAe - QsUs) + Qn-l) 

s=l 

71-1 



(93) 
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as 



E 

p=0 



{A 2 ) p 
AP+ n ~ 2 p\(p + n - 2)! 



p+n-2 



(98) 



1 




n-2 


4n-2 


{dqj _ 




1 




n-2 


4^-2 


fe) _ 





p 1 



oo (A 2 )P 

y 

^ 4Pp\(p + n-2)\ 



e a i 



fi2n-2 



(99) 



so that we finally get the pole 



Pole 



2 2 "- 1 7r 2 e/i- me 

dn 2n-2 ,li=n-l 



n-2 



fi2n - 2 1 n ^ 



(100) 



n g ^+2 (• • • > 
i=i 



which, using yields the contribution to the (3 function in terms of a 2n — 2 dimensional angular 
integral over a finite Euclidean region: 



G 



2m[n 1 ,...,» n ](Pli-)P2m-l) 



(101) 



(n-2)! 



ft, 



n-2 



2n-2 



4rtyj-n+l 
i=n 

Y[G l Vt+2 (...,Ae + q + Q l ^ 1 ) 



j=n-l 

n 

3=1 



i=l 



+ 



cross 



terms 



Expressions (|94D , (|93D and (|101|) for the (3 function of an arbitrary nonlocal coupling are the 
most general results of this paper.! 



9 It is verified easily that with n = 2, to = 2 the (3 function calculation for the nonlocal (j> 4 coupling, as in section 
1, is reproduced. For n = 3, to = 3, and n = l,m = 2, the previous results for the <fi 6 maximally convergent graph 
(2 — 2 — 2) and 4> A tadpole are also verified from the general expression. 
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4 Concluding Remarks 



We have illustrated explicitly how to obtain the renormalization group coefficients in a nonlocal 
scalar field theory. Note that our general results were computed with massive fields, and for this 
reason techniques such as the Gegenbauer polynomial method of |J, valid only for massless propa- 
gators, are not useful here. 

We regard our results as suggestive, but not final. We have obtained integro-differential renor- 
malization group equations by the brute-force technique of exanding, renormalizing and then re- 
summing. We would like to find the rule that allows us to do this in one step, rather than three. 
It would be interesting if, by modifying the analytic structure of the nonlocal theory, we could 
actually give a one-step prescription to isolate the divergent part of the Feynman integral, without 
ever expanding in terms of a formal Taylor expansion. Work along these lines is in progress || . 

We would like to thank Peter Cho, Mike Dugan and Ben Grinstein for numerous discussions. 
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